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Abstract 



A new approach to solving linear ill-posed problems is proposed. 
The approach consists of solving a Cauchy problem for a linear op- 
erator equation and proving that this problem has a global solution 
whose limit at infinity solves the original linear equation. 

1 Introduction 

Let A be a linear, bounded, injective operator on a Hilbert space H, and 
assume that is unbounded and that H^H < \/m, where m > is a con- 
stant. For example, A may be a compact injective linear operator. Consider 
the equation. 



Assume that ( |1 . 1| ) is solvable, so that / = Ay for a unique y & H. Problem 
(|1.1| ) is ill-posed since A~^ is unbounded. Equation ( |LlD cannot be solvable 
for all / G -ff because if A is injective, linear, closed and R{A) = H, then 
A~^ must be bounded (by the Banach theorem). Let fs be given, such that 
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Au = f. 



(1.1) 



\\fs-f\\<S 



(1.2) 



Equation ( |1.1| ) with fs in place of / may have no solution, and if it has a 
solution Us then it may be that \\u — us\\ is large, although J > is small. 
There is a large literature on ill-posed problems since they are important in 
applications. (See e.g.0, 0). In this paper a new approach to solving linear 
ill-posed problems is proposed. This approach consists of the following steps: 

Step 1. Solve the Cauchy problem: 

u = -[Bu + e{t)u- Fs], u{0)=Uo, (1.3) 

where 



u:=—, B:=A*A, Fs:=A*fs, \\Fs - F\\ < 6^^, F = By, 
at 



and 



e{t) G C^[0. oo); e{t) > 0; e{t) \ as t ^ oo; 4^ ^ as t oo. 

(1.4) 

One has ||v4*(/5 — /) || < \/m5., where we have used the estimate \\A\\ = 
\\A*\ \ < vm. 

Examples of functions e{t) satisfying (1.4) can be constructed by the 
formula: ^ 

e{t) = [c+ / h{s)ds]~^, 







where c > is a constant, h{s) > is a continuous function defined for 
all s > 0, such that h{s) — » as s ^ oo and h{s)ds = oo. One has 



\e{t)\ _ 2h{t) 



5 



as t OO. For example, e{t) = i^g^^^2) satisfies (1.4). If 



h{t) 



(2 + t) log(2 + t) 
then 
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This e{t) yields nearly fastest decay of h{t) allowed by the restriction h{s)ds 
oo. 

Step 2. Calculate u{ts), where > is a number which is defined by formula 
([rgp below. 

Then — oo as ^ — and satisfies the inequality: 

\\u{ts)-y\\<v{S) ^0 as 6-^0, (1.5) 

for a certain function r]{6) > 0. If 5 = 0, so that Fs = By, then Step 2 yields 
the relation 

Mm \\u{t) - y\\ = 0. (1.6) 

The foregoing approach is justified in Section 2. Our basic results are 
formulated as follows. 

Theorem 1.1. Assume that equation ( |1.1| ) is uniquely solvable, ( |1.4| ) holds, 
and 6 = 0. Then for any uq, problem (1.3), with F = By replacing Fs, has 
a unique global solution and ( |1.6| ) holds. 

By global solution we mean the solution defined for all t > 0. 

Theorem 1.2. Assume that equation ( |1.1| ) is uniquely solvable, ( |1.4| ) holds, 
and 6 > 0. Then for any Uq problem (|1.3| ) has a unique global solution u{t) 
and there exists a ts ^ oo as 6 ^ 0, such that \\u{ts) — y\\ ^0 as 6 ^ 0. 
The number ts is defined by formula ( |1.9|) . 



;i-7) 



Let y solve (ITT]). Then By = F := A*f and < m. If 

dExy 



m := <PiP,y) ■.= /3 



A + /9 



where Ex is the resolution of the identity of the selfadjoint operator B, 
Ex_Q = Ex, is the minimizer of the function 

h{P,6):=m + ^ (1-8) 
on (0,oo), (see formula (2.20) and Remark 2.3 below), and 
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then ts oo as S ^ 0, rj^S) — > as 5 ^ 0, and 

\im\\uits)-y\\=0. (1.10) 

Because B is injective, zero is not an eigenvalue of B, so, for any y E H, one 
has 1 1 Jq dExvl | — as s — > 0. Therefore ^ as /3 ^ 0, for any fixed 

y. From (2.15) (see below) one gets 

\\u{ts) - y\\ < r]{5) + gsiU) ^ as 5 ^ 0, (1.11) 

where g5{t) is given by the right-hand side of (2.12) with WfsW replacing ||/||. 

Remark 1.1. Theorem 1.2 shows that solving the Cauchy problem (|1.3| ) and 
calculating its solution at a suitable time ts yields a stable solution to ill-posed 
problem (|1.1| ) and this stable approximate solution satisfies the error estimate 

dns). 

For nonlinear ill-posed problems a similar approach is proposed in |T|]. 



2 Proofs 

2.1 Proof of Theorem 1.1 

We start with a simple, known fact: if equation ( |1.1|) is solvable, then it is 
equivalent to the equation 

Bu = A*f = By (2.1) 

Indeed, if Ay = /, then apply A* and get (|2.1|). Conversely, if ( |2.1| ) holds, 
then {B{u — y),u — y) = \\A{u — y)\f = 0, thus Au = Ay and u = y, so ( |1 . 1| ) 
is solvable and its solution is the solution to (|2.1j ). Therefore we will study 
equation ( p.l| ). The operator B = A* A is selfadjoint and nonnegative, that 
is, {Bu,u) > 0. Let Ex be its resolution of the identity. 
We make another observation: If ( |1.4| ) holds, then 

POO 

/ e{t)dt = oo. (2.2) 
Jo 

Indeed, (|1.4]) implies 

e 
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where c = const > 0, so 

d 1 

at e 



e{t) e{0) 
1 



1 



< Co + ct, 



e{t) 

Co := [e{0)]-' > 0, 

and 

sit) > ^ 



Co + ct 

Formula (|2.2| ) follows from the foregoing inequality. 
Consider the problem 

Bw + e{t)w - F = 0, F := A*f = By. (2.3) 

Since B > and e{t) > 0, the solution w{t) of ( |2.3| ) exists, is unique and 
admits the estimate 

M<\\{B + e{t))-'F\\<^^ . (2.4) 
If F = A*f, then (see Remark 2.3 below) one gets: 

\H\ <\\{B + e{t))~-^F\\ = \\{B + e{t))-^A*f\\ < . (2.4') 

2£2 (t) 

Differentiate (p.3|) with respect to t (this is possible by the implicit function 
theorem) and get 

[B + e{t)]w = -iw, ||zi;||<M||^||<^||F||, (2.5) 



where ( p^ ) was used. 
Using (1^) yields: 



kll < — Ikll < (2.5' 
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Denote 



Subtract 



Multiply 



Denote 



z{t) := u{t) - w{t). 
from ( |1.3| ) (with F in place of F^) and get 

z = -w -[B + e{t)]z, z{Q) =UQ-w{<d). 
by z{t) and get 

(i, z) = -{w, z) - {Bz, z) - e{t){z, z). 



urn ■■=9{t) 

Then the inequality {Bz, 2) > and equation ( |2.8| ) imply: 

99 < \\w\\g - e{t)g^. 
Because > 0, it follows from (2.10) and (2.5') that 



a < 



£{t)9{t), 9(0) = \\uo - Wo\ 



so 



9{0) 



2£2 (r) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



Assumption ( |1.4| ) (the last one in (|1.4|)) and ( p.l2| ) imply (use L'Hospital's 
rule) that 

\\u{t) - w{t)\ := git) as t ^ +00. (2.13) 
The existence of the global solution to ( |1.3|) is obvious since equation 



173 ) is linear and the operator B is bounded. 
To prove ( |1.6|) it is sufficient to prove that 

\\w{t) — y|| ^ as t ^ 00. 



(2.14) 
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Indeed, if {^J§ holds then PT^j ) and {^J§ imply: 

\\u{t) - y\\ < \\u{t) - w{t)\\ + \\w{t) - y\\ ^ as t ^ oo. (2.15) 
We now prove (|2.14|) . One has: 

e{t) 



Mt)-y\ 



— — j-dExy- / dExy 
^ + io 



A + £(t) 



Thus 



(2.16) 



(2.17) 



where (l){e,y) := (;/>(£:) is as defined in (|1.7| ). Since B is injective, the point 
A = is not an eigenvalue of B. Therefore 



lim (j){e) = 0, 

e— >0 



(2.18) 



by the Lebesgue dominant convergence theorem. 
Thus ( p.l4| ) follows and Theorem 1.1 is proved. 



□ 



2.2 Proof of Theorem 1.2 

The proof is quite similar to the above, so we indicate only the new points. 
Equation ( p?3|) is now replaced by the equation 



Bw + e{t)w -Fs = 0. 



(2.19) 



Estimates (U), (|2^ , (U), (|23D and (|27[3|) hold with and fs in place 
of F and /, respectively. The main new point is the estimate of w{t) — y: 



\w{t) - y\ 



dExFs 
A + e{t) Jo 
dEx{Fs - F) 



A + £(t) 



< 4>m) + 



dExy 

+mt)) 



(2.20) 



where ||/ — /^H < 5 and estimate (2.4') was used. 
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If (3{6) is the minimizer of the function (1.8), then 

h{p{6),5) := r]{5) ^0 as 5 ^ 0; (3{6) ^ as 6^0. 

The latter relation in (2.21) holds because (p{(3) as /9 — > 0. 
Since e{t) \ as t ^ oo, one can find the unique ts such that 



f2.211 



Thus 



e{ts) = (3{6) ^ as 5 ^ 0. 



\w{ts) - y\\ < r]{6) ^ as 5 ^ 0. 



(2.22) 



(2.23) 

{e, y) does (see 



The function ri{6) = ri{6,y) depends on y because (f){e) = 
formula (|1.7| )). 

Combining ( p.23| ), ( |2.13| ) and ( |2.15| ) one gets the conclusion of Theorem 
1.2. □ 

Remark 2.1. We also give a proof of ( p.l4| ) which does not use the spectral 
theorem. 

From ( p. 3D one gets 



Bx + e{t)x = -e{t)y, x{t) := w{t) - y. (2.24) 
Thus {Bx,x) +e{x,x) = —e{y,x). Since (Bx.x) > and e > 0, one gets 

{x,x) < \{y,x)\, \\x{t)\\ < \\y\\ = const < oo. (2.25) 
Bounded sets in H are weakly compact. Therefore there exists a sequence 



oo such that 



Xfi . x(tn) 



n 



oo 



(2.26) 



where stands for the weak convergence. From ( |2.24D and ( p^.25D it follows 
that 



BXr. 



0, n — »• oo. 



(2.27) 



A monotone hemicontinuous operator is weakly closed. This claim, which 
we prove below, implies that ( |2.26| ) and ( |2.27| ) yield Bx^o = 0. because B is 
injective, x^ = 0, that is, x(t„) ^ 0. From ( 2.25 ) it follows that ||2;(t„)|| 
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as n oo, because {y,x(tn)) — as n — > oo, due to x{tn) 0. By the 
uniqueness of the limit, one concludes that lim = 0, which is (|2.14|) . 

t — ^oo 

Let us now prove the claim. 

We wish to prove that Xn ^ x and Bxn f imply Bx = f provided 
that B is monotone and hemicontinuous . The monotonicity implies {Bxn — 
B{x — ep),Xn — X + ep) > for all e > and all p E H . Take e — > and use 
hemicontinuity of B to get (/ — Bx,p) > ^p E H . Take p = f — Bx to 
obtain Bx = f , as claimed. □ 

The above argument uses standard properties of monotone hemicontinu- 
ous operators j^. 

Remark 2.2. In (|2.23| ) r](5) = 0(6^) is independent of y if y runs through 
a set Sa '■= Sa^R := {y : y = B"-h, \\h\\ < R}, where R > is an arbitrary 
large fixed number and a > 1. //O < a < 1, then ri{6) = 0(6^'^+^), as 6 ^ 0, 
and this estimate is uniform with respect to y E Sa^R. 

Consider, for example, the case a >l.Ify = B^h, then (j){e) in ( ^^201) 
can be chosen for all y G Sa,R simultaneously. Using (|1.7| ), one gets: 

(t){e)=e sup II / dExh\ \ < em''-^R, 

\\h\\<R Jo A + £ 

where a > 1 and e is positive and small. For a fixed 5 > one finds the 
minimizer e{5) = O^S^) of the function -^ + em°-~^R and the minimal value 

ri{S) of this function is 0{S^) . 

If B is compact, then the condition y E Sa means that y belongs to a 
compactum which is the image of a bounded set \ \h\ \ < R under the map B°^. 

The case < a < 1 zs left to the reader. It can be treated by the method 
used above. 

Remark 2.3. It can be checked easily that 

A{A*A + el)-^ = {AA* + eiy^A. 

This implies 

\\iB + eir'A*f\\' = iib + dr%fJ) := J, 
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where B := A* A and h := AA* > 0. 
Thus, 

rm -I 

J = I s{s + e)-'d{ej,f)<-\\f\\^ 

where is the resolution of the identity corresponding to the selfadjoint op- 
erator b. Therefore one gets the following estimate: 

\\{B + eir'A*f\\< ^ 



This estimate was used to obtain estimates (|2.4l ) and ( p. 5' ). For example, 



estimate ( |2.4| ) was replaced by the following one: 



2e2 

and (|2.5|) can be replaced by the estimate: 



B + f/)-'f||<-!T-|l/ll, (2.4') 



IHI<^IIHI<^II/I|. (2.5') 

These estimates were used to improve the estimate for ri{5) in the previous 
remark. 

Estimate (2.4') was used by a suggestion of a referee. The author thanks 
the referee for the suggestion. 

In fact, one can prove a stronger estimate than (2.4'), namely \ \w\ \ < \ \y\\- 
Indeed, multiply (2.3) hj w — y, use the nonnegativity of B and positivity of 
e and get {w,w — y) < 0. Thus < and the desired inequality 

WwW < \ \y\ \ follows. 

Appendix. Let us give an alternative proof of Theorem 1.2. Let us{t) 
solve (1.3), u(t) solve (1.3) with Fs replaced by F, and us{t) and u(t) satisfy 
the same initial condition. Denote ws ■= us{t) — u{t) and let \ \ws\ \ := gsit). 
One has: 

Ws = -[Bws + e{t)ws- hs], 1^5(0) = 0, 

where hs := F^ — F, \\hs\ \ < y/m6 := c6. Multiply the above equation by ws 
in H, use the inequality B > and get 

gs < -^{t)gs + c6. 
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Since ^'^(O) = 0, this implies: 

/■t pt pp ^ 

gs{t) < cS ex.p[— / e{s)ds] / cxp[ / €{s)ds\dp < c—-^. 
Jo Jo Jo 

Thus 

\\us{t)-y\\ < \\us{t) - u{t)\\ + \\u{t)-y\\ < c— + a(t), 

where a{t) :— \ \u{t) — ?/|| — > as t — > oo. Define ts as the minimal minimizer 
of the following function of t for a fixed S > 0: 

c— + a{t) — min :— fJ,{S). 

6[t) 

Since a{t) and e(t) ^ as t ^ oo, one concludes that the minimal 
minimizer oo as 5 — > and fi{6) — >• as (5 ^ 0. Theorem 1.2 is proved 
□ . 
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